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Resumen 


Aldama, A. A. (septiembre-octubre, 2016). Análisis tripartita 
de escurrimientos naturales: aplicación a la caracterización 
de sequías en el río Colorado. Tecnología y Ciencias del Agua, 
7(5), 167-195. 


En este trabajo se presenta el análisis del régimen de 
escurrimiento natural en la cuenca del río Colorado, a la 
luz de lo establecido por el Tratado de 1944 entre México 
y Estados Unidos. Primeramente, se construye el registro 
mensual de escurrimiento natural anual en la estación más 
cercana a la frontera entre ambos países, que corresponde al 
sitio aguas arriba de la presa Imperial. Esto permitió contar 
con valores de septiembre de 1906 a diciembre de 2012, 
de acuerdo con la información más reciente disponible. Se 
calculan los estadísticos más significativos de dicho registro 
y se comprueba que no sigue una distribución normal. Se 
procede a calcular los percentiles 2.28, 15.87, 84.13 y 97.72, 
y con base en ellos se caracteriza el registro histórico en 
periodos extraordinariamente escasos, escasos, normales, 
abundantes y extraordinariamente abundantes. Mediante 
la aplicación de un análisis de Fourier, se determinan los 
periodos dominantes del registro. A través de una novedosa 
descomposición tripartita, se calculan los componentes 
tendenciales, cíclicos dominantes y aleatorios del registro. Se 
demuestra que estos últimos siguen una distribución normal 
y se lleva a cabo la caracterización del registro histórico en 
periodos muy secos, secos, normales, húmedos y muy 
húmedos, de acuerdo con el índice dinámico estandarizado 
de escurrimiento propuesto en este trabajo por el autor. 


Palabras clave: escurrimiento natural anual en la cuenca del 
río Colorado, Tratado de 1944 entre México y Estados Unidos, 
aguas transfronterizas, distribución normal, teorema del 
límite central, percentiles, análisis espectral, descomposición 
tripartita, componente tendencial, componente cíclico, 
componente aleatorio, índice dinámico estandarizado de 
escurrimiento. 


Abstract 


Aldama, A. A. 2016). Three-Part 
Decomposition of Natural Flows: Application to the Colorado River 


(September-October, 


Drought Characterization. Water Technology and Sciences (in 
Spanish), 7(5), 167-195. 


An analysis of the natural flow in the Colorado River basin is 
presented in this paper, within the context of what is established by 
the 1944 Treaty between Mexico and the United States. Firstly, the 
monthly record of natural annual flow at the station closest to the 
border between both countries, which corresponds to the location 
just upstream of the Imperial dam, is built. Thus, values between 
September of 1906 and December of 2012 were included, in 
accordance with the most recently available information. The most 
important statistics of such record are estimated and it is shown 
that the record itself does not follow a normal distribution. The 
2.28, 15.87, 84.13 y 97.72 sample percentiles are computed and 
on the basis of their values, the historical record is characterized in 
terms of extremely scarce, scarce, normal, abundant and extremely 
abundant periods. The dominant periods in the record are identified 
via Fourier analysis. Through a novel three-part decomposition, 
the trend, dominant cyclical and random components of the 
record are identified. It is shown that the last follows a normal 
distribution and a characterization of the historical record in very 
dry, dry, normal, humid and very humid periods is performed, on 
the basis of a dynamic standardized runoff index, proposed by the 
author herein. 


Keywords: Natural annual flow in the Colorado River basin, 
1944 Treaty between México and the United States, transboundary 
waters, normal distribution, Central Limit Theorem, percentiles, 
spectral analysis, three-part decomposition, trend, cyclical 
component, random component, dynamic standardized runoff 
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Introducción 


Los retos que plantea la gestión sustentable 
del agua en las cuencas transfronterizas entre 
México y Estados Unidos son extraordinarios, 
habida cuenta de que el recurso es escaso en 
la región y existe una gran competencia por 
el mismo. Además, la legislación que regula 
el aprovechamiento del líquido en cada país 
y el Tratado sobre Distribución de las Aguas 
Internacionales celebrado entre México y los 
Estados Unidos en 1944 definen un marco 
jurídico complejo que es necesario respetar al 
tomar cualquier decisión relativa al manejo del 
agua. La realidad es que el sistema de cuencas 
transfronterizas, en particular las de los ríos 
Colorado y Bravo, exhibe un frágil equilibrio 
que se trastoca cada vez que ocurre una sequía. 
El Tratado de 1944 menciona ciertas medidas 
especiales de gestión que deben tomarse cuan- 
do ocurra una “extraordinaria sequía”, pero 
desafortunadamente no define el término con 
precisión. Las asignaciones que el Tratado hace 
de las aguas compartidas por las dos naciones 
en las cuencas del Colorado y del Bravo son con- 
dicionantes para el aprovechamiento del recurso 
en ambos países. La interpretación y aplicación 
de las disposiciones contenidas en el "Tratado 
en relación con sequías extraordinarias ha sido 
instrumentada por la Comisión Internacional de 
Límites y Aguas (CILA), de conformidad con las 
atribuciones que le confiere el propio Tratado. 
No obstante, la reciente ocurrencia de periodos 
muy secos tanto en la cuenca del río Colorado 
como en la del río Bravo ha hecho que surjan 
presiones al interior de México y de Estados 
Unidos en el sentido que resulta necesario defi- 
nir un criterio que permita determinar cuándo 
ocurre una extraordinaria sequía y entonces 
aplicar las disposiciones especiales que al res- 
pecto prevé el Tratado. Las secciones mexicana 
y estadounidense de la CILA, con el apoyo de 
las instancias involucradas en la gestión del 
agua en ambos países, entre las que destacan 
la Comisión Nacional del Agua (Conagua) y el 
United States Bureau of Reclamation (USBR), han 
iniciado estudios para caracterizar el régimen 
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de escurrimiento en el río Colorado y contar con 
criterios que permitan, sin ambigiedades, deter- 
minar cuándo es necesario aplicar restricciones 
al consumo del agua en ambos países, a fin de 
que se cumplan las disposiciones generales del 
Tratado. Un avance en este sentido lo constitu- 
yó el Acta 319 (CILA, 2012), firmada el 20 de 
noviembre de 2012, en Coronado, California. A 
través de dicho instrumento jurídico, las seccio- 
nes de la CILA acordaron una serie de medidas 
que incluyen aspectos como la generación y 
conservación de volúmenes adicionales de agua; 
control de salinidad; manejo de variabilidad en 
el abastecimiento de agua en la cuenca; oportu- 
nidad para establecer “agua mexicana creada 
intencionalmente”; intercambio de agua, y 
entrega de agua para el ambiente, entre otras. 
Dicha Acta debe ser reemplazada por un Acta 
Integral a partir del 1? de enero de 2018, lo cual 
subraya la importancia y urgencia de contar con 
los criterios antes referidos para caracterizar el 
régimen de escurrimiento en el río, en particular 
los periodos secos. El propósito de este artículo 
es contribuir al acervo metodológico disponible 
hasta la fecha para llevar a cabo dicha caracteri- 
zación. Para ello, primero es conveniente revisar 
lo dispuesto por el Tratado de 1944. 


Disposiciones relevantes del Tratado de 
1944 


El artículo 10 del Tratado sobre Distribución de 
Aguas Internacionales entre los Estados Unidos 
Mexicanos y los Estados Unidos de América 
establece lo siguiente: 


“De las aguas del río Colorado, cualquiera que 
sea su fuente, se asignan a México: 

a)- Un volumen garantizado de 
1,850,234,000 metros cúbicos (1,500,000 acres 
pies) cada año, que se entregará de acuerdo 
con lo dispuesto en el artículo 15 de este 
Tratado. 

b).- Cualesquier otros volúmenes que 
lleguen a los puntos mexicanos de derivación; 
en la inteligencia de que, cuando a juicio de 


la Sección de los Estados Unidos, en cualquier 
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año exista en el río Colorado agua en exceso 
de la necesaria para abastecer los consumos en 
los Estados Unidos y el volumen garantizado 
anualmente a México de 1,850,234,000 metros 
cúbicos (1,500,000 acres pies), los Estados 
Unidos se obligan a entregar a México, según 
lo establecido en el artículo 15 de este Tratado, 
cantidades adicionales de agua del sistema del 
río Colorado hasta por un volumen total que 
no exceda de 2,096.931,000 metros cúbicos 
(1,700,000 acres pies) anuales. México no 
adquirirá ningún derecho, fuera del que le 
confiere este inciso, por el uso de las aguas 
del sistema del río Colorado para cualquier 
fin, en exceso de 1,850,234,000 metros cúbicos 
(1,500,000 acres pies) anuales. 

En los casos de extraordinaria sequía o de 
serio accidente al sistema de irrigación de los 
Estados Unidos, que haga difícil a éstos entre- 
gar la cantidad garantizada de 1,850,234,000 
metros cúbicos (1,500,000 acres pies), por 
año, el agua asignada a México, según inciso 
a) de este artículo, se reducirá en la misma 
proporción en que se reduzcan los consumos 


en los Estados Unidos”. 


Debe observarse que en el caso de la cuenca 
del río Colorado, el Tratado asigna volúmenes 
de escurrimiento anuales. Por tanto, a fin de 
eliminar el efecto que pueda tener la operación 
de la infraestructura hidráulica construida en 
Estados Unidos, los análisis encaminados a la 
definición del término “extraordinaria sequía” 
deben abordarse en términos de los escurri- 
mientos naturales anuales para el río Colorado. 


El registro mensual de escurrimientos 
naturales anuales en la cuenca del río 
Colorado 


En vista de lo dispuesto por el artículo 10 del 
Tratado, las siguientes secciones estarán enfoca- 
das en analizar el comportamiento estadístico de 
los escurrimientos naturales anuales en la cuenca 
del río Colorado hasta el punto de entrega del 
agua a México por parte de Estados Unidos. 
La Conagua proporcionó al autor el registro de 


escurrimientos naturales mensuales restituidos 
en dicha cuenca, que contiene información de 
octubre de 1905 a diciembre de 2012. Dicho 
registro ha sido integrado por el USBR y es el 
más reciente disponible. 

Sea o, el volumen de escurrimiento natural 
correspondiente al mes m (con m = 1 repre- 
sentando a enero; m = 2, febrero; ..., m = 12, 
diciembre) y al año calendario a, que puede 
tomar valores desde 1905 hasta 2012. Entonces, 
el registro de escurrimientos naturales men- 
suales R, = To, l..:905) + Ú (o... Dome) permite 
construir un registro mensual de escurrimientos 
naturales anuales. Lo anterior, con el propósito 
de capturar el comportamiento detallado de 
la dinámica de escalas de tiempo subanuales 
en los escurrimientos anuales. Para el efecto, 
defínase el siguiente índice anual consecutivo: 


k=a-1905, a4=1905, 1906,..., 2012 (1) 


4 . 1287 
Constrúyase el registro R, =(u,) 2 con ele- 
mentos definidos como sigue: 
(C) = (0) <A 
Uiokem-9 = Ur 2(a-1905)+m-9 =0, uE 1905, 1906,. Mz 2008 


m, a? 


10,11,12 si a=1905 
donde m= 
1,2,...,12 si a=1906,...,2012 (2) 


Entonces, defínase el volumen de escurri- 
miento natural anual correspondiente al índice 
consecutivo j¡ como: 


] 
Usa Y dis] 12 18 pu 1287 (3) 


(=j-11 


El volumen de escurrimiento correspondien- 
te al mes 11 y al año calendario a se obtiene como 


sigue: 

Wa =U racioncrino 51906, 1907 , ss, 2008 
9,10,...,12 si a=1906 

donde m= 


1,2,...,12 si a=1907,...,2012 (4) 


De este modo, se generó el registro men- 
sual de escurrimientos naturales anuales 
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R, =([V,,] 


m,a 1906) ee 


UV. 120007) ,,.1 Nótese que 


1 49=1907 

el registro anterior comienza con el escu- 
rrimiento natural anual del mes de sep- 
tiembre de 1906, pues los volúmenes corres- 
pondientes a los meses anteriores de ese año 
no pueden ser calculados, al no contarse con 
información necesaria. Por tanto, el número 
total de elementos del registro R, es igual a 
12(2012 - 1906) + 4=1 276 =N. En la figura 1 
se muestra la evolución temporal del registro 
mensual de escurrimientos naturales anuales en 
la cuenca del río Colorado. 


Distribución empírica de probabilidad y 
estadísticos principales 


Con el objeto de determinar la distribución 
empírica de frecuencias relativas acumuladas 
del registro mensual de escurrimientos natu- 
rales a numerados consecutivamente, 


Ry= tu, Nes necesario reordenar los elemen- 


0) 


N-1 
cla 


los estadísticos de orden de R,,. En otras pala- 


tos de dicho registro. Al respecto, sean (U, 


bras, (U,,,,)!*, constituye el conjunto de valores 
de los escurrimientos naturales anuales, ordena- 
dos de menor a mayor. pS ahora que el 


k-ésimo valor del conjunto (U o y Corresponde 


(n) 
con el /-ésimo valor de R,,. Entonces, se dice 
que el número de orden o rango de UU, es k + 1. 


Indíquese lo anterior como sigue: 
r(U,)=k+1, si U,=U,,; 
£2=0,1,..,N-1;k=0,1,...,N -1 (5) 


donde r(-) denota el rango. 

Recuérdese ahora que la distribución de 
probabilidad acumulada F(-) de un cierto valor 
X se define como la probabilidad de que X no 
sea excedido por una variable aleatoria muda 
x, esto es: 


F(U) =P(x=<X) (6) 


40000 ; 


35 000 + T 


30 000 + Í 


15000 + 1 | | 


Escurrimiento natural anual (hm?) 


10 000 Í 
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Figura 1. Evolución mensual de los escurrimientos naturales anuales. 
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En congruencia con la ecuación (6), la proba- 
bilidad empírica p, , de que un cierto valor del 
escurrimiento anual U, no sea excedido estará 
dada por: 


_ k+1 HU) 
Pe NA 


(7) 


Con base en la ecuación (7) se obtuvo la 
distribución empírica de probabilidad mostrada 
en la figura 2. 

Es relevante subrayar que en vista de que los 
escurrimientos naturales anuales representan 
la suma de escurrimientos mensuales y éstos 
representan la suma de escurrimientos diarios, 
el teorema del límite central (Meyer, 1973), cuya 
demostración se incluye en Aldama (2011), haría 
pensar que la distribución de aquéllos debería 
aproximarse asintóticamente a una normal. 

Los parámetros estadísticos más importantes 
del registro mensual de escurrimientos naturales 
anuales están incorporados en el cuadro 1. Como 
puede observarse, la mediana es 1.64% inferior 


a la media. Asimismo, el valor del coeficiente de 
variación indica que la variabilidad de los escu- 
rrimientos naturales anuales es moderada. Por 
otra parte, la densidad de probabilidad empírica 
exhibe ligera asimetría positiva (con el máximo 
desplazado hacia la derecha, a diferencia de la 
densidad de probabilidad normal, que es simé- 
trica) y es ligeramente platicúrtica (menos pun- 
tiaguda que la normal). Todo lo anterior puede 
obedecer a efectos asociados con el tamaño de la 
muestra, por lo que no necesariamente desecha 
la posibilidad de que la distribución poblacional 
sea normal, de donde es pertinente evaluar la 
hipótesis de normalidad. 


Prueba de normalidad 


La distribución normal teórica puede ser esti- 
mada como: 


2 


1 u-U 


== 1 u 
AUS 157 3 S 
uy -. u 


du (8) 


Probabilidad empírica 


ad 
0 5 000 10 000 15 000 


20 000 25 000 30 000 35 000 40 000 


Escurrimiento natural anual (hm?) 


Figura 2. Distribución de probabilidad empírica acumulada del registro mensual de los escurrimientos naturales anuales. 
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Cuadro 1. Parámetros más importantes de la distribución de probabilidad empírica del registro mensual 


de escurrimientos naturales anuales. 


Parámetro Valor 
Media 19 971.19 hm? 
Desviación estándar 5 389.66 hm? 
Mediana 19 643.94 hm? 
Coeficiente de variación 0.270 
Coeficiente de asimetría 0.192 
Curtosis -0.559 


donde U y S, representan estimadores de la 
media y la desviación estándar poblacionales 
del registro mensual de los escurrimientos natu- 
rales anuales del río Colorado. Para el cálculo de 
dichos parámetros se han empleado los métodos 
de momentos, mínimos cuadrados y máxima 
verosimilitud. 

En las figuras 3 a 5 se muestran comparacio- 
nes entre las distribuciones empírica y normal 
teórica, estimada por los métodos de momentos, 
mínimos cuadrados y máxima verosimilitud, 
respectivamente. Como puede observarse, los 
datos siguen de manera notable la tendencia de 
la distribución normal. Con el objeto de realizar 
una evaluación cuantitativa, se juzga conve- 
niente aplicar una prueba de normalidad. Entre 
las pruebas de normalidad disponibles, la de 
Kolmogorov-Smirnov (NIST /Sematech, 2006) 
es la menos sensible al tamaño de la muestra, 
por lo que se ha elegido para los efectos de este 
trabajo. 

La hipótesis nula para la aplicación de dicha 
prueba será que la distribución poblacional es 
normal. A partir de los datos de la muestra se 
construye el estadístico siguiente: 


D = sup E (U;N)-F, (U¡U, S,) (9) 
u 
donde F (UN) representa la extensión continua 


de la distribución de probabilidad empírica 
acumulada, definida como: 


1 N 
A (10) 
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donde Lise es la función indicadora, que toma 
el valor de 1 si U,,,< U y 0 si U,,,> U. Obsérvese 
que F (UN) está relacionada con la probabili- 
dad empírica mediante la siguiente expresión: 


E[U,. ¿Nel =p,;k=1,2,..., No (11) 


(e) 

Se procedió a calcular el parámetro D y se 
obtuvo el siguiente resultado: D = 0.0464, con 
parámetros estimados mediante el método de 
momentos D = 0.0445, con parámetros estima- 
dos mediante el método de mínimos cuadra- 
dos, y D = 0.0465, con parámetros estimados 
mediante el método de máxima verosimilitud. 
El valor crítico para este parámetro, corres- 
pondiente a un nivel de significancia de 5%, es 
D,, =1.36/,[N. =0.0381 (West Chester Univer- 
sity, 2009). Dado que en los tres casos, D > DD, 
se descarta la hipótesis nula. En conclusión, 
los datos analizados no siguen una distribución 
normal. Este resultado llama poderosamente 
la atención, en vista de la expectativa creada 
por la aplicación del teorema del límite cen- 
tral, mismo que predice que datos como los 
analizados se deberían distribuir en forma 
asintóticamente normal. Sobre este particular 
se abundará después. 


Caracterización del régimen mensual de 
escurrimiento natural anual por medio 
de percentiles 


Sea x una variable aleatoria normalmente dis- 
tribuida, con media u y desviación estándar o, 
de modo que su distribución de probabilidad 
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Probabilidad 


0.0 5 000.0 10000.0 15000.0 20000.0 25000.0 30000.0  35000.0  40000.0 


Escurrimiento natural anual (hm?) 


Figura 3. Comparación entre las distribuciones empírica y normal teórica, estimada por el método de momentos, del registro 
mensual de escurrimientos naturales anuales. 


E) «==» Normal teórica (mínimos cuadrados) a 
WM 


Probabilidad 


0.0 5000.0  10000.0 15000.0  20000.0 25000.0 30000.0  35000.0  40000.0 


Escurrimiento natural anual (hm?) 


Figura 4. Comparación entre las distribuciones empírica y normal teórica, estimada por el método de mínimos cuadrados, 
del registro mensual de escurrimientos naturales. 


Tecnología y Ciencias del Agua, vol. VIL, núm. S, septiembre-octubre de 2016, pp. 167-195 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VI, núm. 5, septiembre-octubre de 2016, pp. 167-195 


Aldama, Análisis tripartita de escurrimientos naturales: aplicación a la caracterización de sequías en el río Colorado 


«=== Normal teórica (máxima verosimilitud) 


+ Empírica 


Probabilidad 
o 
u 


0.0 5000.0  10000.0  15000.0 


20000.0 25000.0  30000.0 35000.0  40000.0 


Escurrimiento natural anual (hm?) 


Figura 5. Comparación entre las distribuciones empírica y normal teórica, estimada por el método de máxima verosimilitud, 


del registro mensual de escurrimientos naturales. 


esté dada por: 


5x0) fexp 2) de (12) 


Introdúzcase la variable estandarizada z, 
definida como: 


gane (13) 


Sustituyendo la ecuación (13) en la ecuación 
(12) se obtiene la así denominada distribución 
normal estándar: 


0:0.)= fexp 36 ja (14) 


Por definición, la media de z es nula y la 
variancia unitaria. 
Nótese que, de acuerdo con la ecuación (13): 
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z=-2=>x=u-20 
z=-1>x=u-0 
z=1>x=u+0 
z=2=>x=u+20 (15) 
Dado que la distribución normal estándar 
está tabulada, se conocen los valores de probabi- 


lidad que corresponden a distintos valores de la 
variable estandarizada z. En efecto, se sabe que: 


z=-2 > Q(-2;0,1)= P(z <-2) = 2.28% 
2=-1> Q(-1;0,1)= P(z <-1)=15.87% 


z=1> 0(1;0,1)=P(z <1) = 84.13% 
2=2=>0(2;0,1)=P(z<2)=97.72% (16) 
De acuerdo con la ecuación (4.5), z = -2 


representa el percentil 2.28; z = -1, el percentil 
15.87; z= 1, el percentil 84.13; y z =2, el percentil 
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97.72 de la distribución normal estándar. Con 
base en las ecuaciones (13) y (16) es posible 
escribir que para una variable aleatoria x que 
esté normalmente distribuida: 


P(x <u-20)=P(z < -2) = 2.28% 
P(u-20<x=<u-0)=P(-2<z<-1)=13.59% 
P(u-o<x<u+0)=P(-1<z <1)=68.26% 
Plu+o<x<u+20)=P(l<z=<2)=13.59% 
P(x>u+20)=P(z>2)=2.28% (17) 


Las probabilidades indicadas en la ecuación 
(17) se ilustran gráficamente en las figuras 6 a 
10. Debe subrayarse que las interpretaciones 
expresadas por la ecuación (17) y en las leyen- 
das de las figuras 6 a 10 sólo son válidas para 
variables normalmente distribuidas. 

No obstante, el concepto de percentil es ge- 
neral y se puede manejar en forma no paramé- 
trica, esto es, sin que esté referido a distribución 
teórica alguna. En efecto, sea y una variable 


aleatoria con distribución de probabilidad cual- 
quiera W(y) = P(n < y), donde y es una variable 
muda. Entonces, por definición, el percentil 
«€ [0,1] (en fracción) es el valor y, tal que W(y,) 
= P(n < y) = k. En consecuencia, en general la 
probabilidad de que se presente un valor menor 
o igual al percentil 2.28 es 2.28%; se presente un 
valor entre los percentiles 2.28 y 15.87 es 13.59%; 
se presente un valor entre los percentiles 15.87 
y 84.13 es 68.26%; se presente un valor entre los 
percentiles 84.13 y 97.72 es 13.59%, y se presente 
un valor mayor al percentil 97.7 es 2.28%. 

Con base en la anterior discusión, tiene 
sentido definir un evento aleatorio con valor: 
“extraordinariamente bajo” al que sea menor 
que el percentil 2.28; valor “ordinariamente ba- 
jo” al que esté entre los percentiles 2.28 y 15.87; 
“normal” al que esté entre los percentiles 15.87 y 
84.13, “ordinariamente alto” al que esté entre los 
percentiles 84.13 y 97.72; y “extraordinariamente 
alto” al que sea mayor que el percentil 97.72. Pa- 
ra pronta referencia, a los percentiles 2.28, 15.87, 
84.13 y 97.72 se les denominará percentiles índice. 


Densidad de probabilidad 


0 1 2 3 4 lo 


Variable estandarizada 


Figura 6. P(x < u - 20) = P(z<-2). 
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Densidad de probabilidad 


Variable estandarizada 


Figura 7. P(u-20<x<u-0)=P(-2<z<-1). 


Densidad de probabilidad 


Variable estandarizada 


Figura 8. P(u-o <x<u+0)=P(1<z<1). 
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Figura 9. P(u+0 <x<u+20)=P(l<z<2). 
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Variable estandarizada 


Figura 10. P(x > u +20) =P(z > 2). 
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Los percentiles índice empíricos de la distri- 
bución de probabilidad del registro mensual de 
escurrimientos anuales mostrada en la figura 2 
se tienen en el cuadro 2. En el cuadro 3 se pro- 
pone una clasificación del registro mensual de 
escurrimientos anuales basada en los percentiles 
índice contenidos en el cuadro 2. 

En la figura 11 se muestra una representación 
gráfica de la clasificación de los escurrimientos 
naturales anuales registrados en el periodo 
de septiembre de 1906 a diciembre de 2012, 
utilizando el sistema de clasificación de escurri- 
mientos presentado en el cuadro 3. De acuerdo 
con la clasificación propuesta en el cuadro 3, 
dicho sistema, con base en la información dispo- 
nible en el registro empleado, estadísticamente 
2.35% del registro mensual de escurrimientos 
naturales anuales, se clasificaría como “extraor- 
dinariamente escaso”; 13.56%, “escaso”; 68.18%, 
“normal”; 13.56%, “abundante”, y 2.35%, “ex- 
traordinariamente abundante”. Las pequeñas 
diferencias entre estos porcentajes observados y 
los teóricos implicados en el cuadro 2 obedecen 
a efectos de redondeo. 


Determinación de los ciclos naturales 
del régimen mensual de escurrimientos 
naturales anuales 


De conformidad con la notación introducida 
en la ecuación (3), sea R,, =(U;) el registro 
mensual de escurrimientos naturales anuales 
identificados mediante el índice ¡ que adopta 
valores enteros progresivos de 0 a N - 1. Tam- 
bién recuérdese que dicho indicador se relaciona 


con el mes 71 y el año a del siguiente modo: 


n=12(4-1906)+m-9,4=1906,1907,..., 2012 


9,10,...,12 si a=1906 
donde m= 
1,2,...,12 si a=1907,...,2012 (18) 


La autocovariancia del registro mensual 
de escurrimientos naturales anuales se estima 


como sigue: 
1 N-m-1 =Ñ —N 
Ruumn => — u-4u yu, -u ) 
uu ,m N -m 2 ( j ) ( jam ) 
m=0,1,..,M (19) 


Cuadro 2. Percentiles índice empíricos de la distribución de probabilidad del registro mensual de escurrimientos naturales 


Percentil Valor (hm?) 
2.28 9 218.8 
15.87 14 394.8 
84.13 26 196.3 
97.72 30 848.3 


Cuadro 3. Clasificación de escurrimientos naturales anuales en el río Colorado. 


Clasificación Porcentaje de probabilidad Intervalo de escurrimiento anual (hm') 
Extraordinariamente escaso 2.28 Menor o igual que 9 218.8 
Escaso 13.59 Mayor que 9 218.8 y menor o igual que 14 394.8 
Normal 68.26 Mayor que 14 394.8 y menor o igual que 26 196.3 
Abundante 13.59 Mayor que 26 196.3 y menor o igual que 30 848.3 
Extraordinariamente abundante 2.28 Mayor que 30 848.3 
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Figura 11. Clasificación de los escurrimientos naturales anuales registrados en el río Colorado. 


donde: 


u, (20) 


es la media muestral de los escurrimientos natu- 
rales anuales, y M << N es un apropiadamente 
elegido valor máximo del sesgo discreto (por lo 
general como una fracción de N). 

Por otra parte, la autocovariancia atempera- 
da se estima como: 


Ruum="0,R yy ¿m=0,L..,M (21) 
donde 1,, es el m-ésimo valor de la ventana de 
Hann, dada por la siguiente expresión: 


Za cos SY m=0,1,..,P 
W = 


m 


0 ¿m=P+1,..,N-1 (22) 


La selección de la ventana de Hann para 
estimar la autocovariancia atemperada se fun- 
damenta en el análisis realizado por Aldama 


(2010), quien demostró que es la que produce 
una estimación más suave o estadísticamente 
estable del espectro, sin dificultar la clara iden- 
tificación de las frecuencias dominantes. 

El espectro discreto se calcula como la 
transformada de Fourier discreta de la auto- 
covariancia atemperada dada por la ecuación 
(21). En la práctica se emplea el algoritmo de 
la “transformada rápida de Fourier” (Stearns, 
1975), mismo que funciona con series cuyo nú- 
mero total de datos es igual a una potencia de 2. 
Por lo anterior, se acostumbra agregar ceros al 
final de la serie de valores de la autocovariancia, 
hasta alcanzar la potencia de 2 más cercana. 
Además, por definición, la autocovariancia es 
una función par, lo cual se traduce en que el 
espectro sea real y par. En términos del cálculo 
de la transformada discreta de Fourier, esto es 
equivalente a repetir los valores de la imagen 
especular de la autocovariancia. En términos 
precisos, se emplea la autocovariancia atempe- 
rada extendida, estimada como sigue: 
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Ruun ¿m=0,1,..., M 


Ruum=10 ;m=M+1,M+2,...M' M'=2% R= pa 


In2 


Ruuzr=n;m=M"+1,M'+2,..,2M'-1 (23) 


donde el símbolo [-] se denomina “techo” y se 
interpreta como el número entero más pequeño 
que no es menor que su argumento. 

Así, el espectro discreto se estima como 


sigue: 


an 2M'-1_, . 
Suun= y Ruune mM; y =0,1,...,2M'-1 (24) 


m=0 


donde At = 1 mes, dado que se trabaja con el 
registro mensual de escurrimientos naturales 
anuales. Además, en vista de que el espectro 
es par y la transformada discreta de Fourier es 
periódica por construcción, sólo son relevantes 


los valores para n =0, 1,..., M, que corresponden 
a los valores originales de la autocovariancia 
atemperada (sin ceros adicionados). Los valores 
(positivos) respectivos de frecuencia y periodo 


son: 
2nn TN 
o, = = ¿n=0,L..,M (25) 
2M'M MA 
a LM (26) 
n 


Se aplicaron las ecuaciones (19), (21) y (23) 
para calcular la autocovariancia original, la au- 
tocovariancia atemperada y la autocovariancia 
atemperada extendida usando M= 600 (= N / 2), 
por lo que M*=1 024 y 2M*= 2 048. En la figura 
12 se muestra la autocovariancia atemperada del 
registro mensual de los escurrimientos naturales 
anuales. 
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Figura 12. Autocovariancia del registro mensual de escurrimientos anuales del río Colorado, atemperada 


con la ventana de Hann. 


ISSN 0187-8336 


Aldama, Análisis tripartita de escurrimientos naturales: aplicación a la caracterización de sequías en el río Colorado 


El espectro asociado con el registro mensual 
de escurrimientos naturales anuales se estimó 
haciendo uso de la ecuación (24) y su represen- 
tación gráfica se muestra en la figura 13. Como 
se puede observar ahí, se identifican claramente 
tres rangos en el comportamiento del espectro: 
(1) de altas amplitudes; (2) autosimilar o fractal, 
en el que el espectro sigue una ley de potencia 
negativa, y (3) disipativo. En la figura 14 se real- 
zan con círculos trazados con línea discontinua 
los máximos locales asociados con las frecuen- 
cias dominantes en el espectro, mismas que se 
identifican de modo cuantitativo en el cuadro 4, 
junto con los periodos que les corresponden. En 
dicho cuadro se muestran periodos dominantes 
de alrededor de 60 años y de aproximadamente 
7a 13 años, y de 2a 5 años. Dichos periodos se 
han ilustrado en la figura 15, sobre la gráfica de 
los escurrimientos naturales anuales. 

Como se puede observar, el periodo de 60 
años se asocia con ciclos largos o súperdece- 
nales, mismos que de manera aproximada se 


asocian con la recurrencia de eventos extraor- 
dinariamente escasos o extraordinariamente 
abundantes. El periodo de 7 a 13 años se relacio- 
na con ciclos de duración intermedia o decenal, 
que pueden o no cambiar la clasificación del 
régimen de escurrimiento, y el de 3 a 5 años 
se vincula con ciclos cortos o subdecenales que 
por lo general no cambian la clasificación del 
régimen de escurrimiento (esto es, todo el ciclo 
se mantiene como normal, abundante o escaso). 


Un nuevo criterio dinámico para 

la caracterización del régimen de 
escurrimiento de un río: aplicación al río 
Colorado. 


En vista de que los escurrimientos naturales 
anuales en la cuenca del río Colorado repre- 
sentan la suma de escurrimientos mensuales 
y éstos, a su vez, la suma de escurrimientos 
diarios, así como lo establecido por el teorema 
del límite central, sería de esperarse que dichos 
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Figura 13. Espectro correspondiente al registro mensual de escurrimientos naturales anuales en la cuenca del río Colorado. 
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Cuadro 4. Frecuencias y periodos dominantes del registro mensual de escurrimientos naturales anuales. 
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Frecuencia (1/mes) Periodo (meses) Periodo (años) 
0.0092 682.7 56.7 
0.0399 157.5 13.1 
0.0736 85.3 7.1 
0.1012 62.1 5.1 
0.1166 53.9 4.5 
0.1963 32.0 2.7 
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Figura 14. Identificación de frecuencias dominantes en el registro mensual de escurrimientos naturales anuales 


en la cuenca del río Colorado. 


escurrimientos anuales siguieran una distribu- 
ción normal. No obstante, se ha demostrado que 
ello no es así. En esta sección se busca entender 
el porqué de este hecho y se propone una me- 
todología alterna de análisis del régimen de los 
escurrimientos naturales anuales. 

Antes de proceder, es pertinente destacar 
que el método de clasificación de escurrimientos 
en términos de percentiles, mismo que es más o 
menos convencional y que ha sido presentado 
antes en este trabajo, corresponde a un criterio 
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estático debido a que utiliza estadísticos que 
aplican por igual, con independencia de la 
fecha en la que haya ocurrido un valor deter- 
minado de escurrimiento natural anual que se 
desee clasificar. De hecho, podría decirse que 
prácticamente todos los criterios existentes para 
clasificar el régimen de variables climáticas o 
hidrológicas son estáticos (como los criterios 
incluidos en el monitor de sequías de América 
del Norte; véase por ejemplo, Aldama, 2013). 
Esto no toma en cuenta que el clima de una re- 
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Figura 15. Periodos subdecenales, decenales y superdecenales identificados en el registro mensual de escurrimientos naturales 


anuales en la cuenca del río Colorado. 


gión y el régimen hidrológico de una cuenca son 
fenómenos dinámicos que poseen variabilidad 
de pequeña, mediana y gran escala. Al respecto, 
resulta conveniente recordar las definiciones 
que propone el Glosario Hidrológico Interna- 
cional (OMM/UNESCO 2012) para “sequía”, 
como “ausencia prolongada o escasez acusada 
de precipitación”; para “sequía hidrológica”, 
como “periodo de clima anormalmente seco 
lo suficientemente prolongado para ocasionar 
una disminución apreciable en el caudal de 
los ríos, nivel de los lagos y/o un agotamiento 
de la humedad del suelo y un descenso en los 
niveles de aguas subterráneas por debajo de 
sus valores normales”, y para “clima”, como 
“síntesis de las condiciones meteorológicas 
en un lugar determinado, caracterizada por 
estadísticas a largo plazo (valores medios, 
varianzas, probabilidades de valores extremos, 
etc.) de los elementos meteorológicos en dicho 
lugar”. Las definiciones anteriores ponen de 


manifiesto que en el estudio del régimen hidro- 
lógico de un río es necesario tener en cuenta las 
estadísticas de largo plazo (algo que sí hacen 
los criterios estáticos), la desviación respecto 
de “valores normales” (ídem), y la duración de 
los diferentes estados que puede exhibir dicho 
régimen (lo cual, por definición, no pueden 
hacer los criterios estáticos). Queda claro que 
en el análisis hace falta involucrar a un prota- 
gonista fundamental: el tiempo. Al respecto, la 
experiencia ha demostrado que, además de que 
pueden existir efectos de largo plazo (lo que ha 
dado en llamarse “cambio climático”), el clima 
y la hidrología de una cuenca tienden a exhibir 
ciertos patrones cíclicos, determinados por la 
dinámica diurna y estacional del movimiento 
de la tierra con respecto al motor de los fenó- 
menos climático-oceánico-hidrológicos, esto es, 
el sol, así como la dinámica anual y multianual 
de ese astro. A lo anterior se pueden sumar va- 
riaciones de pequeña escala que podrían estar 
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determinadas por causas que no son fácilmente 
discernibles y que sería conveniente caracterizar 
como inciertas, aleatorias o fortuitas. 
Tomando en cuenta la discusión anterior, 
defínase V, de la siguiente forma: 
V, = Ma (27) 
donde, como se recordará, V, representa el 
escurrimiento natural anual correspondiente al 


mes 1 del año calendario a, y A representa el 
año fraccionario, dado por: 


A=a-1+m/12 (28) 


La ecuación (28) implica que A coincide con a 
en el mes de diciembre de cada año calendario. 

Considérese la siguiente descomposición 
tripartita del escurrimiento natural anual co- 
rrespondiente año fraccionario A: 


V, =V +V¿+V/ (29) 


donde V; representa la tendencia lineal del es- 
currimiento natural anual; V;, el componente 
cíclico, y v/, el componente aleatorio o fortuito 
de dicho escurrimiento. 

A diferencia del análisis estadístico tradi- 
cional (i.e., criterios estáticos), que se basa en 
suponer que todo el escurrimiento natural 
anual es una variable aleatoria, en la ecuación 
(29) se postula que el escurrimiento está com- 
puesto por una componente determinista dada por 
V, + Ví y una componente aleatoria o fortuita, 
dada por V/. Con fundamento en el teorema 
del límite central, si la ecuación (29) es correcta, 
la componente aleatoria V/ poseerá una distri- 
bución normal. Lo anterior, en vista de que la 
descomposición (29) es lineal y homogénea, por 
lo que es de esperarse que tanto la tendencia 
como el componente cíclico y el componente 
aleatorio puedan expresarse como la suma de 
sus correspondientes contrapartes para el es- 
currimiento natural diario para todo el año. A 
continuación se presenta una metodología para 
determinar cada una de las tres componentes 
del escurrimiento natural anual que intervienen 
en la ecuación (29). 
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La ecuación (29) es el modelo más simple o 
parsimonioso (en el sentido del principio de la 
parsimonia o “Navaja de Ockham”, Audi, 1999) 
que captura las características deseables en un 
criterio de caracterización dinámico, ya que, 
como quedará claro en la líneas subsecuentes, 
el elemento tiempo interviene de manera directa 
en la determinación tanto de la componente 
tendencial como la cíclica. 

La tendencia del escurrimiento natural 
anual (V/) se obtiene por simple regresión, 
como se acostumbra hacerlo en análisis de se- 
ries de tiempo (Falk et al., 2006). Utilizando los 
datos del registro de escurrimientos naturales 
Ry = MV 212-1906) a U rado ) ió se obtiene 
la siguiente expresión para la tendencia del 
escurrimiento natural anual: 


Vi=0A4+p (30) 


donde a. = -41.6166697876027 hm*/año y 
p = 101 493.046478986 hm. En la figura 16 se 
muestra una representación gráfica de la ten- 
dencia de escurrimiento natural anual en el río 
Colorado. 

Ahora bien, defínase el primer residuo de 
escurrimiento natural anual, va, como la dife- 
rencia entre el escurrimiento natural anual y su 
tendencia, esto es: 


VO=V, -V;=V¡+V/ (31) 


donde en la última igualdad se ha empleado la 
ecuación (29). En analogía con las ecuaciones 
(3) y (4), teniendo en mente la ecuación (28), se 
pueden definir los elementos del conjunto de 
primeros residuos identificados con un índice 


: (1) N-1 : . 
consecutivo, tU; ) ¡o , COMO sigue: 


¡0 yu 


12(a-1906)+m-9 = Y a m/12/ 


a=1906,1907,..., 201: 


9,10,..., 12 si a=1906 
donde m= 
1,2,...,12 si a=1907,...,2012 (32) 


Exprésense dichos elementos mediante una 
representación discreta de Fourier, esto es: 
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Figura 16. Tendencia del registro mensual de escurrimientos naturales anuales en la cuenca del río Colorado. 


N-1 
up” - NUPEINYA;j=0,L..., N-1 (33) 
k=0 


donde los coeficientes de Fourier ((U% E ,) están 


dados por la transformada discreta de Fourier 
Mpal A 
j Hi? 


de los primeros residuos ((U 
mo mo N-1 . . 
qe E zu Mad = y IEA 


=0 


bill Nel (34) 


donde la tilde y E [-] denotan la transforma- 
da discreta de Fourier y Af es el intervalo de 
tiempo que separa a los distintos valores de 
los escurrimientos naturales anuales (en el 
presente caso, un mes). La frecuencia a la que 
corresponde el valor de U'” es: 


En 180) 
NAt T, 
donde Aw representa el intervalo de frecuencia 


que separa los distintos valores de U0, y T, el 


periodo asociado con W Evidentemente, dichos 
parámetros están dados por las siguientes ex- 


presiones: 
qa (36) 
NAf 
7, k=0,L,., N-1 (87) 


Las ecuaciones (33) y (34) constituyen el par 
discreto de Fourier de los primeros residuos de 


los escurrimientos naturales anuales. 
Al 
0 


estimador-base de su espectro (Stearns, 1975), 


El periodograma de (Uj”)2;' constituye el 


y se define como: 
a 40 = UP IÚPT;k=0,L..,N-1 (38) 
donde el asterisco denota conjugado complejo. 


En este trabajo se postula que el componente 
cíclico del escurrimiento natural anual (V;) 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VII, núm. 5, septiembre-octubre de 2016, pp. 167-195 


Tecnología y Ciencias del Agua, vol. VI, núm. 5, septiembre-octubre de 2016, pp. 167-195 


Aldama, Análisis tripartita de escurrimientos naturales: aplicación a la caracterización de sequías en el río Colorado 


corresponde a la reconstrucción de Fourier que 
contenga los modos dominantes del periodo- 
grama 9. 

Los elementos del conjunto de componentes 
cíclicos identificados con un índice consecutivo 
(u; 14) se definen como sigue: 


Us 


12(a-1906)+m-9 


=V" 


a-1+m/12 / 


a=1906,1907,..., 200 


9,10,...,12 si a=1906 
donde m= 
1,2,...,12 si a=1907,...,2012 (39) 


De manera similar, los elementos del 
conjunto de componentes aleatorios iden- 
“ps » . . N-1 
tificados con un índice consecutivo (qu ) j=0 ) se 


definen como sigue: 


Sj a , 
9,10,...,12 si 4a=1906 

donde m= 

1,2,...,12 si a=1907,...,2012 (40) 


uf 


12(a-1906 )+m-9 


a=1906,1907,..., 2008 


Ahora bien, los modos dominantes de Fou- 
rier son aquéllos que corresponden a los valores 
mayores de %,. A continuación se desarrolla un 
método para determinarlos. El valor máximo 
del periodograma está dado por: 


Prngx = MÉX IZ, (41) 
Defínase el periodograma condicional (2, %"“) 
de la siguiente manera: 

RO ALA PVP (42) 
donde 0 < p< 1 es la fracción de Y”, , que 
define la condición de corte de las amplitudes 
del periodograma, y .47[:] es la función escalón 
unitario de Heaviside, dada por: 


Osi T<0 


A El 


(43) 


Como puede entonces observarse, la 
ecuación (42) deja intactos los valores del pe- 
riodograma Y, que son mayores o iguales que 
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el valor de corte p%,, , haciendo cero el resto. 
Evidentemente, si se permitiere a p asumir el 
valor de cero, el periodograma condicional 
sería igual al periodograma original. De hecho, 
cuanto más cercana esté p de cero, tanto más 
parecido será 2% ah. 

Supóngase que el valor de p fuere conocido, 
entonces el componente cíclico del escurri- 
miento natural anual (u)) estaría dado por la 
reconstrucción condicional de Fourier, que invo- 
lucra sólo a los nodos supervivientes del tamiz 
involucrado en la definición del periodograma 
condicional a través de la ecuación (42), esto es: 


N-1 


UA EIA 
k=0 
=D) Li N=1 (44) 


Por las razones antes expuestas, mientras 
más pequeña sea p, más cercano será u a un" 
y, por tanto, más pequeño será U/, pues de 
conformidad con la ecuaciones (31) y (32): 


Y = ur Un (45) 


La descomposición tripartita del escurri- 
miento anual, dada por la ecuación (29), prescribe 
que dicha variable está dada por la suma de un 
componente determinista, dado por la suma de la 
tendencia y el componente cíclico, y un compo- 
nente aleatorio. Por tanto, si el primero crece, el 
segundo debe menguar. Pero es precisamente el 
componente aleatorio ul el que permite carac- 
terizar probabilistamente la incertidumbre en el 
comportamiento de los escurrimientos naturales 
anuales. A fin de no sesgar la descripción del 
régimen de escurrimiento a través del modelo 
representado por la ecuación (29), se buscará 
maximizar la variancia explicada por el com- 
ponente aleatorio, ul . No obstante, en vista de 
lo establecido por el teorema del límite central, 
dicho componente debe poseer una distribución 
normal. 

Tomando en cuenta lo anterior, se propone 
la determinación de los componentes cíclico y 
aleatorio del primer residuo de escurrimientos 
naturales anuales a través de la solución del 
siguiente problema de optimización no lineal: 
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p'=máxp 
O<p<l 


tal que: 
Nc 


U¡= ) AR, JUDEA; j=0,L,..., N=1; y 


máx 
k=0 


u/=u-U;;j=0,1....,N-1 posea una distribución 
normal (46) 


Se aplicó el método descrito al análisis del 
registro mensual de escurrimientos naturales 
anuales del río Colorado. En la figura 17 se 
muestra una representación gráfica de los 
primeros residuos, calculados por medio de la 
ecuación (31). En otras palabras, dicha figura 
indica el registro mensual de escurrimientos 
naturales anuales sin tendencia (por tanto, su 
media es nula). 

Es pertinente hacer notar que para la apli- 
cación práctica de las ecuaciones (33) a (46) es 
necesario incrementar el número de elementos 


del conjunto (U$”)*;', de modo que sea posible 
el uso de la transformada rápida de Fourier. Es- 
to implica la inclusión de elementos nulos de la 
siguiente forma: ur =0paraj=N,N +1,...,N*, 
donde N* es la potencia de 2 más cercana (desde 
arriba) a N. En este caso, entonces, se sustituye, 
donde sea aplicable, N =1 276 por N*=2 048 en 
las ecuaciones (33) a (46). 

Se resolvió numéricamente el problema de 
optimación planteado por la ecuación (46) y se 
obtuvo que el valor óptimo de la fracción de cor- 
te para el registro mensual disponible es p =0.4. 
Se calculó el componente cíclico mediante la 
reconstrucción condicional de Fourier expresa- 
da por la ecuación (44), y el segundo residuo 
o componente aleatorio, mediante la ecuación 
(45). Debe hacerse notar que las frecuencias 
que contienen el componente cíclico del primer 
residuo son esencialmente las mismas que se 
determinaron mediante el análisis espectral 
realizado antes. 
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Figura 17. Primer residuo de escurrimiento natural anual en el río Colorado. 
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En la figura 18 se muestra una comparación 
entre el primer residuo y el componente cíclico. 
Puede observarse que el componente cíclico 
tiene un comportamiento similar al del primer 
residuo, pero que sus oscilaciones son más 
suaves y de menor amplitud. Por otra parte, en 
la figura 19 se muestra el componente aleatorio. 
Se observa que la variabilidad de éste es signi- 
ficativamente menor que la del primer residuo. 
De hecho, hay una disminución de 49.5% en la 
variancia. 

Resta demostrar que el componente aleatorio 
se distribuya de manera normal. En las figuras 
20 a 22 se muestran las comparaciones entre la 
distribución empírica de probabilidad de dicho 
componente y las distribuciones normales teó- 
ricas, con parámetros estimados mediante los 
métodos de momentos, mínimos cuadrados y 
máxima verosimilitud. Como puede observarse, 
el ajuste es muy satisfactorio en los tres casos. 

En términos cuantitativos, se aplicó la 
prueba de normalidad de Kolmogorov- 


Smirnov, empleando la ecuación (9), y se ob- 
tuvieron los siguientes resultados: D = 0.0365, 
con parámetros estimados mediante el méto- 
do de momentos; D = 0.0370, con parámetros 
estimados mediante el método de mínimos 
cuadrados; D = 0.0364, con parámetros es- 
timados mediante el método de máxima 
verosimilitud. Dado que en los tres casos, 
D < D,, = 0.0388, se concluye que con un nivel 
de significancia de 5%, el componente aleato- 
rio se distribuye de forma normal. Esto con- 
firma la expectativa planteada al inicio, a la 
luz de lo establecido por el teorema del límite 
central. En vista de que el ajuste es mejor, en 
términos de D, para los parámetros estimados 
mediante el método de máxima verosimilitud 
(además de que arroja estimadores consisten- 
tes; Hogg, McKean, € Craig, 2012), se decidió 
emplear los valores correspondientes para la 
media (U/) y desviación estándar (S,y), mis- 
mos que resultaron ser O hm* (como era de 
esperarse) y 3 731.18 hm”, respectivamente. 
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Figura 18. Comparación entre el primer residuo y el componente cíclico del escurrimiento natural anual del río Colorado. 
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Figura 19. Componente aleatorio del escurrimiento natural anual del río Colorado. 
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Figura 20. Comparación entre la distribución empírica de probabilidad del componente aleatorio y la distribución normal con 
parámetros estimados empleando el método de momentos. 
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Se propone ahora un nuevo índice para clasi- 
ficar hidrológicamente el régimen de escurrimien- 
tos naturales de una cuenca. La denominación 
respectiva es la de “índice dinámico estandari- 
zado de escurrimiento” o DSRI (abreviatura de 
Dynamic Standarized Runoff Index, en analogía 
al SPI, abreviatura de Standarized Precipitation 
Index, y de SDI, Standarized Drought Index; véase 
Nalbantis € Tsakiris, 2009). 

El DSRI se define como la diferencia entre el 
componente aleatorio del escurrimiento natural 
y su media, normalizada al dividirse entre la 
desviación estándar correspondiente. Emplean- 
do la notación de los componentes aleatorios del 
escurrimiento natural anual en el río Colorado, 
se tiene que: 

-U y 


uf 
DSRI/? = a (47) 


uf 


En el cuadro 5 se muestra la clasificación 
propuesta en términos del DSRI. 


Debe recordarse que el componente aleatorio 
es la diferencia entre el escurrimiento natural 
anual y la suma de la tendencia y el componen- 
te cíclico. Por tanto, DSRI = k significa que el 
escurrimiento difiere en k veces la desviación 
estándar del componente aleatorio de dicha su- 
ma. En la figura 23 se muestra una comparación 
entre el escurrimiento natural y la suma de la 
tendencia, y el componente cíclico, misma que 
corresponde a DSRI=0. 

El DSRI reconoce que la condición de nor- 
malidad no se puede caracterizar con un valor, 
sino con un rango. Por otra parte, dado que la 
suma de la tendencia y el componente cíclico 
es variable en el tiempo, también lo serán los 
límites entre los diversos rubros de clasificación 
contenidos en el cuadro 7. Es esto lo que da el 
carácter dinámico al DSRI. Nalbantis y Tsakiris 
(2009), así como otros autores, han asignado 
probabilidades a rangos de índices empleados 
por ellos con base en la distribución normal sin 
demostrar que los datos que han empleado se 
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Figura 21. Comparación entre la distribución empírica de probabilidad del componente aleatorio y la distribución normal 


con parámetros estimados empleando el método de mínimos cuadrados. 
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Figura 22. Comparación entre la distribución empírica de probabilidad del componente aleatorio y la distribución normal 
con parámetros estimados empleando el método de mínimos cuadrados. 


Cuadro 5. Clasificación de escurrimientos en términos del índice dinámico estandarizado de escurrimiento. 


Criterio Clasificación Probabilidad (%) 
DSRI>2 Muy húmedo 2.28 ES 
1<DSRIs 2 Húmedo 13.59 S 
-1< DSRI<1 Normal 68.26 5 
-2<DSRIs-1 Seco 13.59 S 
a 
DSRI < -2 Muy seco 2.28 2 
distribuyan de modo normal. A diferencia de 1944. Puede observarse que la condición de F 
lo anterior, el DSRI garantiza por definición normalidad corresponde a la zona entre DSRI = E 
que los datos a los que se le aplica sean nor- -1 y DSRI= 1, lo cual reconoce que el clima y sus E 
malmente distribuidos. Las probabilidades efectos tienen variaciones periódicas, además de a 
consignadas en el cuadro 5 son congruentes con ciclos de gran escala que pueden manifestarse El 
esta observación, dado que corresponden a las como tendencias. 5 
probabilidades de la componente aleatoria del De acuerdo con el DSRI, en el periodo de E 
escurrimiento natural anual. octubre de 1906 a diciembre de 2012 se presentó a 
En las figura 24 se muestra la clasificación de 3.292% de meses muy secos, 10.266% de meses 3 
escurrimientos naturales anuales en la cuenca secos, 67.555% de meses normales, 17.633% de 
del río Colorado, hasta el punto de entrega de meses húmedos y 1.254% de meses muy hú- 3 
agua a México, en los términos del Tratado de medos. Debe aclararse que estos porcentajes ya S 
e 
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Figura 23. Comparación entre el registro mensual de escurrimiento anual en el río Colorado y la suma de su tendencia y 


componente cíclico. 
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Figura 24. Clasificación de escurrimientos naturales anuales en la cuenca del río Colorado. 
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incluyen el efecto de la tendencia. Los periodos 
muy secos se presentaron de agosto de 1931 a 
febrero de 1932, marzo de 1935, mayo de 1935, 
de julio a diciembre de 1963, de julio de 1977 a 
abril de 1978, de julio de 1981 a marzo de 1982, 
y de agosto de 1994 a febrero de 1995. 

El método propuesto en este trabajo, basado 
en el DSRI, constituye un sistema de clasifica- 
ción objetivo del régimen de escurrimiento 
natural en un río y es consistente con la distri- 
bución normal de probabilidad, incorporando 
asimismo el concepto de normalidad variable, 
que reconoce que el clima y sus efectos exhiben 
comportamiento periódico y variaciones a largo 
plazo. Por otra parte, la metodología de análisis 
desarrollada en este capítulo reconoce que en 
el registro de escurrimientos naturales pueden 
estar involucrados componentes no aleatorios, lo 
cual implicaría que los mismos no sean inde- 
pendientes, hecho que explica el porqué si los 
valores de los escurrimientos no se descompo- 
nen en componentes deterministas y compo- 
nentes aleatorias no se satisfaga la normalidad 
distribucional, en vista de que entonces se viola 
una de las hipótesis fundamentales del teorema 
del límite central, que es la de independencia 
estadística de los valores muestrales. 


Comparación de los criterios estático y 
dinámico para la clasificación de sequía 


Resulta interesante comparar el criterio estático, 
basado en percentiles, y el criterio dinámico, 
sustentado en la descomposición tripartita, que 
han sido presentados en las secciones anteriores 
para caracterizar el régimen mensual de escu- 
rrimiento natural anual. Dicha comparación se 
presenta en el cuadro 6. Puede observarse que 


el criterio del DSRI implica que los meses muy 
secos son un poco más frecuentes; los secos, 
un poco más infrecuentes; los normales, de 
frecuencia similar; los húmedos, un poco más 
frecuentes, y los muy húmedos son un poco 
más infrecuentes. Esto se debe en particular 
a la influencia de la tendencia y componente 
cíclico, que caracterizan la dinámica propia de 
la variabilidad del escurrimiento natural anual 
en la cuenca del río Colorado. 

Por la importancia que tienen en el manejo 
sustentable del agua en la región bajo estudio, 
resultan de particular interés los periodos 
escasos y extraordinariamente escasos. En la 
figura 25 se muestra una comparación de los 
criterios basados en percentiles y el DSRI para 
determinar la ocurrencia de dichos periodos. 


Conclusiones 


El Tratado sobre Distribución de las Aguas 
Internacionales celebrado entre México y Esta- 
dos Unidos en 1944 establece la obligación que 
tiene el segundo país de asignar al primero un 
volumen anual de 1 850.234 hm? de las aguas 
del río Colorado. Asimismo, implica que cuando 
ocurra una extraordinaria sequía en la parte de 
la cuenca ubicada en los Estados Unidos, que 
haga difícil a éstos entregar dicha cantidad a 
México, el volumen entregado se reducirá en 
la misma proporción en que disminuyan los 
consumos en los Estados Unidos. Desafortuna- 
damente, el Tratado no contiene una definición 
de “sequía extraordinaria”. Los gobiernos de 
México y Estados Unidos, en el seno de la CILA 
han iniciado una colaboración para caracterizar 
el régimen de escurrimiento en el río. El autor 
de este artículo ha asesorado técnicamente a la 


Cuadro 6. Comparación de criterios para la clasificación del régimen mensual de escurrimiento anual en la cuenca del río 


Colorado. 
EE % de meses muy % de meses % de meses % de meses % de meses muy 
Criterio A 4 
secos secos normales húmedos húmedos 
Percentiles 2.351 13.558 68.182 13.558 2.351 
DSRI 3.292 10.266 67.555 17.633 1.254 
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Figura 25. Comparación de criterios para la caracterización de periodos secos y muy secos. 
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Conagua y a la Sección Mexicana de la CILA 
en dicha materia. En vista de que el Tratado 
habla de volúmenes anuales, se ha abordado la 
caracterización del régimen de escurrimientos 
naturales anuales, a fin de eliminar los efectos 
propios de la operación de infraestructura 
hidráulica. En este trabajo se presenta un aná- 
lisis estadístico detallado de dicho régimen. En 
particular, se demuestra que los datos crudos 
del régimen mensual de escurrimiento natural 
anual no siguen una distribución normal, 
resultado que llama la atención en vista de 
que el teorema del límite central predice que 
una variable aleatoria, definida como la suma 
de un número significativo de otras variables 
aleatorias, tiende asintóticamente a distribuirse 
de forma normal. También se presenta el análisis 
espectral de los datos y se identifican periodos 
dominantes de oscilación del orden de 60 años, 
que representan ciclos superdecenales, entre 7 y 
13 años, que representan ciclos decenales; y de 
2 a 5 años, que significan ciclos subdecenales. 
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También se propone un criterio estático de cla- 
sificación del régimen de escurrimiento, basado 
en percentiles. La propuesta más significativa de 
este trabajo consiste en presentar el desarrollo 
de una nueva metodología para clasificar regí- 
menes de escurrimiento en general y sequías 
en particular. Dicha metodología se basa en la 
descomposición tripartita del escurrimiento 
natural anual en una componente tendencial, 
una componente cíclica y una componente 
aleatoria. Se presenta el procedimiento para 
determinar el componente cíclico en términos 
de un periodograma condicional, basado en la 
determinación de los principales modos de Fou- 
rier, sujeto a que el componente aleatorio posea 
una distribución normal. Asimismo, se explica 
la razón por la cual los escurrimientos crudos 
no siguen dicha distribución, pues contienen 
componentes que no son estadísticamente 
independientes, condición sine qua non para 
la aplicación del teorema del límite central. Se 
propone, asimismo, el índice dinámico estanda- 
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rizado de escurrimiento (DSRI, por sus siglas en 
inglés), que provee un medio objetivo, riguroso 
y racional para clasificar escurrimientos, en 
tanto que está basado en el teorema del límite 
central e incorpora los efectos dinámicos ten- 
denciales y cíclicos de un registro hidrológico. El 
uso combinado y juicioso del DSRI y del criterio 
estático basado en percentiles puede contribuir 
a la gestión sustentable del agua transfronteriza 
en la cuenca del río Colorado, permitiendo en 
especial hacer previsiones para los periodos 
secos y muy secos. 
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